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Abstract 

The energy level statistics of 2D electrons with spin-orbit scattering are 
considered near the disorder induced metal-insulator transition. Using the 
Ando model, the nearest-level-spacing distribution is calculated numerically 
at the critical point. It is shown that the critical spacing distribution is size 
independent and has a Poisson-like decay at large spacings as distinct from 
the Gaussian asymptotic form obtained by the random-matrix theory. 
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Typeset using REVTeK 



Recently, the statistical description of electronic properties at the critical region near the 
Anderson transition in disordered systems has been the subject of great interest. Universal 
behavior was proposed for the variance of the ensemble averaged number of energy 
levels in a given interval, {{6Ny) = {N"^) — (N)"^, and for the nearest-level-spacing distribu- 
tion P{s). A relation {{6N)'^) = a (N) could be extracted from numerical calculations for 
the three-dimensional Anderson model, while for P{s) a universal combination of the well 
known Wigner surmise (small spacings) and the Poisson distribution (large spacings) was 
proposed. Recent analytical calculations suggested, however, another universal relation 
{{6Ny) = a(ip{N)"' / (3 near the critical transition (mobility edge) and the corresponding P{s) 
was different from both the Poisson and the Wigner limit. Universal means here that the 
result depends only on the symmetry class /3, the correlation exponent v = [d{l —7)]"^, and 
the spatial dimension d. The statistics of energy levels are governed by the symmetry of the 
Hamiltonian belonging to certain universality classes of the corresponding Gaussian ensem- 
bles of random matrices: orthogonal (/3 = 1), unitary (/3 = 2) and symplectic (/5 = 4) 

In two-dimensional (2D) disordered systems a complete Anderson transition is expected 
only for the symplectic symmetry class. Two different models for describing the localization 
problem in two dimensions in the presence of spin-orbit coupling were proposed For 
the Ando model |0] the critical behavior of the localization length was previously studied by 
means of the transfer-matrix method The statistics of energy levels for symplectic 



symmetry in connection with the results of the random-matrix theory (RMT) was considered 



in the Evangelou-Ziman model [|T^] and the distribution P{s) at the transition point for this 



model has already been mentioned ||13[ . The multifractal properties of the critical eigenstates 



in the symplectic case were investigated in |]l4| , p!5| . 

In this paper we study the level statistics of electrons with spin-orbit coupling in 2D 
disordered systems. By using the Ando model we calculate numerically the correlations in 
the exact one-electron energy spectrum at the metal-insulator transition. Our main result 
is that the probability density P{s) of neighboring levels in the critical region is universal, 
i.e. it does not depend on the size of the system, and has a novel form. 



It is known that in the metallic phase P{s) is very close to the Wigner surmise |TG 



appropriate for the symplectic symmetry class of the random Hamiltonians |T7 



Pgse{s) = ^s^ exp {-^s^) , (1) 

where s is measured in units of the mean level spacing. In the strongly localized regime the 
energy spectrum is completely uncorrelated, and the distribution of the level spacings obeys 
the Poisson law 

Pp(s) = exp(-s). (2) 

In addition to these two universal distributions one can expect that there is a third form 
of P{s), which corresponds exactly to the metal- insulator transition. This new universal 
statistics was already found for the 3D disordered systems without spin-orbit scattering 
(/3 = 1) and confirmed numerically PJT8[]. An analogous result was recently obtained in the 



unitary case {(3 = 2) |T^, where time-reversal symmetry is broken by the magnetic field. 

It is reasonable to suppose that a similar scale-invariant universality of P{s) also holds at 
the critical point for 2D disordered electrons with spin-orbit coupling. Thus, we expect that 
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in the thermodynamic hmit there exist three possible hmiting situations for (3 = 4 namely, 
the Wigner surmise Eq. ([l|) for the metallic regime, the Poisson law Eq. (j^) for the insulating 
regime, and the critical distribution P{s) at mobility edge. Therefore, if L — oo, the level 
statistics change discontinuously twice, from the delocalized regime to the transition point 
and, then, from the transition point to the localized regime. 

In order to calculate the critical level statistics we start with the Hamiltonian of the 
Ando model 

+ y{n,a]m,a')cl^c^^,, (3) 

^,cr n,m,cr,a' 

where c]^^ and c^^ are the creation and annihilation operators of an electron at a lattice site 
n = (x, y) with the spin component a, and m denotes the sites adjacent to the site n. The 
on-site energy e„ is randomly distributed around zero according to a box distribution with 
a width W. The parameter W specifies the degree of the disorder. The transfer matrices 
V{n, a; m, a') = Vx,Vy depend on the direction 



V _ ^1 '^2 V - 

V x,y;x+l,y T/ T/ ' -t/ \^) 





and describe the hopping between the nearest neighbor sites in the lattice. The strength 
of the spin-orbit coupling is given by the parameter S = V2/V, where V = {V^ + Vg^)^*^^ 
is taken as the unit of energy. In what follows we consider 5 = 1/2. It was earlier found 



by the transfer-matrix method |jTO| that the metal-insulator transition in the middle of the 
band e/V = occurs at a disorder Wc/V = 5.74. 

Applying periodic boundary conditions, the exact discrete eigenvalue spectrum has 
been obtained from a numerical diagonalization of the Hamiltonian Eq. @ using a Lanc- 
zos algorithm. For square lattices of linear size L/a = 50 (100) energy intervals [—1,0] 
([—0.5,0]) were chosen. Within these intervals the density of states is almost constant, 
p = (AL^)^^ ^ 0.127. Here, a is the lattice constant and A the mean level spacing. The 
eigenvalues taken from these intervals belong to the critical region where the correlation 
length ^ is larger than the system size L. The total numbers of eigenvalues N were 94672 
(300 realizations) and 101744 (160 realizations) for L = 50 a and L = 100 a, respectively. 

Fig. m displays the level spacing distribution function P{s) calculated for the different 
sizes of the system at the metal-insulator transition. One can see that independent of 
the system size the data lie on a common curve. This critical P{s) is very close to Pgse 
and passes the point Sq ~ 1.63, where the two limiting distributions -Pgse(s) and Pp{s) 
cross. For small spacings P{s) oc s^. It is interesting to notice that our results differ from 
those obtained previously for smaller systems of the Evangelou-Ziman spin-orbit coupling 
model p[. 

We applied the fitting function which was recently proposed for the description of the 
critical level distribution [Q, 

P{s) = Bs^ exp{-As^-^) (5) 

and found A = 2.77 ± 0.05, B = 17.8 ± 0.8 and 7 = 0.28 ± 0.03. Using a confidence level 
of 95% (a = 0.05) the fitting procedure within the range < s < 3 where the statistics of 
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our numerical data is rather good yielded = 29.2 which is less than the expected value 
Xa = 63.4. Hence, the analytical formula (Eq. |^) can be accepted within an approximate 
relative error of ix^/NY^^ ^ 2%. Although, in the range < s < 3 the calculated Pfs) is 
in good agreement with Eq. (^, the obtained exponent 7 gives, however, a different value 
of the correlation length exponent u = [(1 — 7) rf]^^ ^ 0.7 as compared to the numerically 
obtained value u = 2.75 using the transfer-matrix method PJTO|]. 

Moreover, for large spacings our results deviate markedly from the above Eq. (inset 
of Fig. p. Instead, the behavior of P{s) in the range of the spacings 1.5 < s ^ 4 is well 
described by the Poisson-like asymptotic form 

P{s) oc exp(— Ks), (6) 

with a coefficient k = 4.0 ± 0.2. This decay is much slower than both the Gaussian decay 
of Eq. (0) and the intermediate decay from Eq. (|^), but faster than that for the insulating 
regime Eq. (Q). A similar exponential tail of the critical P{s) was found also in the 3D-case 
without spin-orbit interactions [^] for which the coefficient n is approximately a factor of 2 
less. Such an asymptotic behavior of the critical P{s) is in good agreement with a suggestion 
about the Gaussian form for the distribution of the number of levels lying within a given 
energy interval p|,pl|. 

In order to diminish the magnitude of relative fluctuations due to the limited number 
of realizations and to analyze the asymptotic behavior of the level spacing distribution in 
detail, it is more convenient to consider a total probability function I{s) = P{s') ds'. 
This quantity implies a portion of spacings which are larger than a given s. It is clear that 
/(O) = 1 and /q°° ds' = 1 regardless of the disorder. In the strongly localized regime 
/p(s) = exp(— s), and Igse{s) can be calculated from the RMT. The results of the numerical 
calculations for the critical I{s) are shown in Fig. ^ One observes again a discrepancy from 
both the GSE asymptotic and the /(s) obtained from Eq. (|^), particularly when s is large. 
The form of the critical /(s) is not sensitive to the change of the lattice size, in analogy with 
P{s). The Poisson tail as described by Eq. (H) is recovered for large spacings: ln/(s) oc —s. 
We have checked that the slope of the linear behavior of ln/(s) at the transition does not 
depend on the width of the energy interval from which the levels are taken, as long as they 
belong to the critical region. In the range of very small probability (s > 3) the larger 
fluctuations observed are due to insufficient statistical data. However, the accuracy of the 
total calculated I{s) is higher than that of P{s). One should notice that the behavior of the 
obtained level spacing distribution does obviously not reflect the information on the critical 
exponent u and the dimensionality d in the form as it is expressed by Eq. (^. 

In conclusion, we have presented results of computer simulations of the nearest-level- 
spacing distribution P{s) for 2D disordered systems in the presence of spin-orbit interaction. 
Exactly at the metal-insulator transition the P{s) and, consequently, the total probability 
of neighboring spacings I{s) do not depend on the system size and are different from the 
universal limiting distributions corresponding to the metallic and the insulating regime. 
They appear to exhibit critical behavior at the disorder Wc/V = 5.74 and finite-size scaling 
properties around the critical point. The obtained large-s parts of P{s) and I{s) are shown 
to have a Poisson-like decay so that InP(s) = —n s where k ~ 4.0 is larger as compared to 
the insulating regime {k = 1). 
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FIGURES 



FIG. 1. The level spacing distribution P{s) for two different system sizes at the critical disorder 
Wc/V = 5.74. Solid lines correspond to PcsEis) and Pp{s), respectively. The dotted line is a fit 
of Eq. (^). The inset shows the large-s behavior of P{s). The straight line is a fit according to 
Eq. (I). 

FIG. 2. The integrated probability I{s) at the critical disorder Wc = 5.74. The solid curves 
are IgseI-s) and /p(s) for the metallic and insulating phases, respectively. The dotted line is 
obtained from Eq. (^). The inset shows the large-s part of I{s). The straight line fitting the data 
is ln/(s) = -4.0 s + 3.6. 
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